Draft version February 5, 2008 

Preprint typeset using I£T^X style cmulatcapj v. 9/08/03 



EVOLUTION OF X-RAY CLUSTER SCALING RELATIONS IN SIMULATIONS WITH RADIATIVE COOLING 

AND NON-GRAVITATIONAL HEATING 

Orrarujee Muanwong 1 , Scott T. Kay 2-3 and Peter A. Thomas 3 

Draft version February 5, 2008 

ABSTRACT 

We investigate the redshift dependence of X-ray cluster scaling relations drawn from three hydrody- 
namic simulations of the ACDM cosmology: a Radiative model that incorporates radiative cooling of 
the gas, a Preheating model that additionally heats the gas uniformly at high redshift, and a Feedback 
model that self-consistently heats cold gas in proportion to its local star-formation rate. While all 
three models are capable of reproducing the observed local Lx-Tx relation, they predict substantially 
different results at high redshift (to z = 1.5), with the Radiative, Preheating and Feedback models 
predicting strongly positive, mildly positive and mildly negative evolution, respectively. 
The physical explanation for these differences lies in the structure of the intracluster medium. All 
three models predict significant temperature fluctuations at any given radius due to the presence of 
cool subclumps and, in the case of the Feedback simulation, reheated gas. The mean gas temperature 
lies above th e dynamical temperature of the halo for all models at z = 0, but differs between models 
at higher redshift with the Radiative model having the lowest mean gaswos temperature at z = 1.5. 
We have not attempted to model the scaling relations in a manner that mimics the observational 
selection effects, nor has a consistent observational picture yet emerged. Nevertheless, evolution of 
the scaling relations promises to be a powerful probe of the physics of entropy generation in clusters. 
First indications are that early, widespread heating is favored over an extended period of heating that 
is associated with galaxy formation. 

Subject headings: galaxies: clusters: general, cosmology: theory 



1. INTRODUCTION 

X-ray scaling relations of galaxy clusters, namely the 
temperature-mass, Tx-M, relation and the luminosity- 
temperature, Lx-Tx, relation, play a pivotal role when 
using the ab undance of clusters to constrain cosmological 
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Allen et al.ll2003t iHenrv II2004D . It is well known, how- 



ever, that accurate calibration of scaling relations is cru- 
cial to avoid a major source of systematic error. For 
example, the Tx—M relation is widely used by many 
of these authors to constrain the amplitude of mass fluc- 
tuations, conventionally defined using the parameter, a$. 
Systematic deviations in the normalization of the Tx—M 
relation, partic ularly due to how c luster mass is esti- 
mated (e.g. see iHorner et al. lll99 9) is amplified by the 
steep slope of the te mperature fun ction, leading to large 
variations in as (see IHenrv 1120041 for a discussion of re- 
cent results). 

As far as the Lx-Tx relation is concerned, the dis- 
crepancies are more prominent as Lx is highly sensi- 
tive to the thermodynamics of the of the inner intra- 
cluster medium (ICM), and can yield different values for 
both normalizations and slopes ijEdge^^Sjewartl [19911 
IWhite. Jones fc Formanl Il997t kllen fc Fabian I Il99st 

1 Department of Physics, Faculty of Science, Khon Kaen Uni- 
versity, Khon Kaen, 40002, Thailand; orrmua@kku.ac.th. 

2 Astrophysics, Denys Wilkinson Building, Keble Road, Univer- 
sity of Oxford, Oxford, OX1 3RH, UK. 

3 Astronomy Centre, Department of Physics and Astronomy, 
School of Science and Technology, University of Sussex, Falmer, 
Brighton, BN19QH, UK. 



IMarkevitch 11799% 1x7 



Wu 1 120001) . The situation is 
further complicated by the fact that clusters do not scale 
sclf-similarly, as would be the case (approximately) if 
the only sour ce of heating was via gravitational infall 
l)Kaiser I IT986). This makes the problem more difficult 
to investigate theoretically, although it allows studies of 
cluster scaling relations to reveal more information on 
the physics governing the structure of the intracluster 
medium. 

The departure from self-similarity can be attributed 
to an increase in the entropy of the gas that particu- 
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formed to investigate the effects of various physical 
processes that can raise the e ntropy of the gas, based 
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Measurements of how cluster scaling relations evolve 
with redshift allow even tighter constraints to be placed 
on cosmological parameters (and entropy generation 
models), and observations of cluster properties at high 
redshift are now starting to become available, owing pri- 
marily to the high sensitivity of Chandra and XMM- 
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Newton. From a theoretical point of view, this is an 
exciting phase as we can now fully exploit the availabil- 
ity of our simulated distant clusters and compare their 
X-ray properties with real observations. It is therefore 
timely to investigate further the effects of entropy gener- 
ation on the evolution of clu ster scaling relations as the 
available data for high-redshift systems accumulates. 

In this paper, we will use cosmological hydrodynam - 
ical simulations described in Muanwo ng et aLl (2002), 
hereafter MTKP02, and in lKav et al. I <|2004() . hereafter 
KTJP04, to trace the evolution of the cluster popula- 
tion to high redshift (z — 1.5). Our results will primar- 
ily focus on three (Radiative, Preheating and Feedback) 
models, all able to reproduce the local Lx-Tx relation. 
The aims of this paper are to determine how the scaling 
relations evolve with redshift in the three models and 
to discover what the evolution of scaling relations can 
teach us about non-gravitational processes occurring in 
clusters. 

The rest of this paper is outlined as follows. In 
Section [21 we introduce the X-ray scaling relations and 
summarize our present observational knowledge of these 
quantities. Details of our simulated cluster populations 
are presented in Section |3| In Section 0] we present our 
main results, first at z — 0, where the models are in good 
agreement with each other and the observations, then as 
a function of redshift, where the models predict widely 
different results. We discuss the implications of these 
differences in Section and demonstrate that the degree 
of X-ray evolution is driven by the supply of cold, low 
entropy gas. Finally, we summarize our conclusions in 
Section |6] 



2. X-RAY CLUSTER SCALING RELATIONS 

Kaiser (1986) derived the following relations for tem- 
perature 

T x cx Afi (1 + z), (1) 

and luminosity 



L x ctM3 (1 + z)5 
ocT x 2 (l + z)i 



(2) 
(3) 

assuming the distribution of gas and dark matter in clus- 
ters is perfectly self-similar and the X-ray emission is 
primarily thermal bremsstrahlung radiation. Observed 
clusters do not form a self-similar population but it is 
nevertheless convenient to describe their behavior using 
a generalized power-law form 



Y = C (z)X a = Y X a (l + z) 



(4) 



where Cq(z) and Y determine the normalization, a is 
the slope of the relation (in log-space) and A determines 
how the relation evolves with redshift. Our main results 
will focus on the determination of A. 

Observationally, attempts to measure the Tx— M rela- 
tion at high redshift are currently in their infancy, as they 
require temperature profiles to be measured so that their 
mass can be estimated, but initia l results are consistent 
with self-similar evolution (A ~ l. lMaughan et al~1l2005t 
iKotov fc Vikhlinin 1120051) . 

Measuring the Lx-Tx relation at higher redshift is a 
somewhat sim pler prospect, and has been attempted by 
many authors ijMushotzkv fe Scharf l!99^:lFairlev et al. I 
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Fig. 1. — Evolution in the Lx-Tx relation as measured from 
various high redshift cluster samples. 
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We summarize recent results that adopt a low-density 
flat cosmology in Figure attempting to include in the 
size of the error bars the uncertainty in A due to the 
choice of local relation (when quoted by the authors). 
Although the present situation is by no means clear, 
taking all results at face value generally favors positive 
evolution (0 < A < 2) with the latest results being 
consistent with sel f-similar evolution (A = 3/2). Larger 
samples of high redshift clusters (such as that expected 
from the XMM-Newton Cluster Survey, IRomer et al. I 
120011) will be crucial to accurately constrain the degree 
of evolution in the Lx-Tx relation. 

3. SIMULATED CLUSTER POPULATIONS 

Our results are drawn from three similarly-sized N- 
body/SPH simulations of the ACDM cosmology, which 
have already been published in MTKP02 and KTJP04. 
The simulation box in MTKP02 has a comoving side of 
100/i _1 Mpc with 160 3 particles each of gas and dark 
matter, whose particle masses are set to 2.6 x 10 9 and 
2.1 x 10 10 h^ 1 M , respectively. The box used in KTJP04 
is bigger with a side of 120/i _1 Mpc using 256 3 particles 
each of gas and dark matter, whose particle masses are 
1.3 x 10 9 and 7.3 x 10 9 h' 1 M , respectively. Full details 
can be found in the articles. The key difference between 
the simulations is the model used to raise the entropy of 
the intracluster gas, summarized as follows: 

1. A Radiative model where the excess entropy origi- 
nated from the removal of low entropy gas to form 
stars, causing higher entropy gas to flow adiabati- 
cally into the core from larger radii (MTKP02). 
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2. A Preheating model where entropy was generated 
impulsively by uniformly heating the gas by 1.5 keV 
per particle at z = 4 (MTKP02). 

3. A Feedback model where the entropy of (on aver- 
age) 10 per cent of cooled gas in high density re- 
gions was raised by 1000 keV cm 2 , mimicking the 
effects of heating due to stars and active galactic 
nuclei (KTJP04). 

These three models differ in the timing and distribu- 
tion of entropy generation in the intracluster medium. 
The Radiative model has no explicit feedback of energy 
but relies on the removal of low-entropy gas via cool- 
ing; as such it represents a minimal heating model. The 
Preheating model contains distributed heating at high 
redshift such as might occur if entropy generation occurs 
mainly in low-mass galaxies. By contrast heating in the 
Feedback model occurs solely in high-density regions. In 
all our models, there is very little star formation before a 
redshift of z=4 after which the star-formation rate (sfr) 
begins to rise rapidly. In the Preheating simulation the 
sfr is then strongly suppressed, whereas in the other two 
simulations it peaks at a redshift of z=2 and then de- 
clines back down to low values by the present day with 
a time- variation that matches that of the star-formation 
history of the Universe. The global baryon fraction in 
stars (and cold gas) at z=0 is 0.002, 0.076 and 0.127 in 
the Preheating, Radiative and Feedback simulations, re- 
spectively. The largest of these corresponds to a stellar 
mass density of fi*=0.006; thus none of the models has 
excessive star-formation. These models are far from ex- 
haustive and their precise details should not be taken too 
seriously. The purpose of this paper is not to examine 
particular models but to illustrate that the evolution of 
the X-ray scaling relations can provide a powerful dis- 
criminant between different classes of model. 

3.1. Cluster identification and properties 

Clusters were selected at four rcdshifts (z =0, 0.5, 1 
& 1.5) using the procedure outlined in MTKP02. They 
are defined to be spheres of matter, centered on the dark 
matter density maximum, with total mass 



properties averaged within -Rsoo- The mass, 



M A = -ttR 3 a A Pc0 {1 + z) 3 



(5) 



where p c o = 3Hq /8ttG is the critical density at z = 0. We 
set A = 500 as it corresponds to a sufficiently large ra- 
dius such that the results are not dominated by the core, 
as well as corresponding approximately to the extent of 
cu rrent X-ray observations. Furth ermore, as was shown 
bv lRowlev. Thomas fc KaTI |20?)1 . the X -ray properties 
of simulated clusters within an overdensity of 500 exhibit 
less scatter than within the virial radius. Our choice of 
scaling with redshift 4 is independent of cosmology and 
would allow the simple power-law scalings to be recov- 
ered (equations Ill2l fel3l if the clusters were structurally 
self-similar. 

We consider scaling relations involving mass, three 
measures of temperature, and luminosity, for particle 

4 Many authors prefer to adopt the redshift scaling of the critical 
density, E(z) 2 = n m (l + z) 3 + S^a (for a flat universe), rather than 
the background density, (1 + z) s . 
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E 
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(6) 



where the sum runs over all particles, of mass m,. The 
dynamical temperature, 



kT c 



dyn 



(7) 



where a = (2/3)/i?7iH ~ 4. 2x 10~ 16 keV for a fully ionized 
primordial plasma, assuming the ratio of specific heats 
for a monatomic ideal gas, 7 = 5/3, and the mean atomic 
weight of a zero metalicity gas, /imn = 10~ 24 g. The 
first sum in the numerator runs over all gas particles, 
of temperature, Ti, whereas the second sum runs over 
particles of all types, of speed v, as measured in the center 
of momentum frame of the cluster. 

We also consider the mass- weighted temperature of hot 
(T > 10 5 K) gas, 



kTp- 



hot ' 



i.kTi 



Si, hot m i 



(8) 



and we approximate the X-ray temperature of a cluster 
using the bolometric emission-weighted temperature, 



kTho\ 



Si,hot m l p l k h rt(T i , Z)Tj 

Ei,hot TO ift A bol(?i,Z) 



(9) 



where p, is the density and Aboi is the bolometric cooling 
funct ion used in our simulations ({Sutherland fc Dopita I 
1993); for the Radiative and Preheating runs, Z = 
0.3(t/t )Z Q (MTKP02), and for the Feedback run, Z = 
0.3Z Q (KTJP04). Finally, the X-ray luminosity is ap- 
proximated by the bolometric emission-weighted lumi- 
nosity 

■m l p l K ho[ (T i , Z) ^ 

i,hot 



^bol 



(/xm H ) 2 



It has been shown recently that the emission-weighted 
temperature is not an accurate diagnostic of cluster tem- 
perature, overpredicting the spectroscopic temperature 
by ~ 20 — 30 per cent when th e emission is predom- 
inantly thermal b remsstrahlung ( Mazzo tta et al. 1 120041: 
IRasia et al. II2005J) . At lower temperatures (kT < 3keV), 
line emission from heavy elem ents makes the p roblem sig- 
nificantly more complicated (Vikhlinin 2005). The vol- 
ume sampled by our simulations (~ 100 h Mpc) means 
that we have very few clusters with T > 3keV, and so a 
more accurate measure of the cluster temperature would 
require significantly more effort than applying a simple 
formula to our data. We therefore leave such improve- 
ments to future work, when larger samples of simulated 
clusters are available. It would not affect the conclusions 
of this paper. 

3.2. Cluster catalogues 

Tabled lists the numbers of clusters in our catalogues 
for each of the simulations at all 4 redshifts. The first 
row for each model gives the total number of clusters 
in our catalogues, down to a minimum mass, M500 = 
1.2 x 10 13 h~ x M©, corresponding to ~ 500 dark mat- 
ter particles in the Radiative and Preheating simulations, 
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Table 1 

Numbers of clusters at various redshifts 



Redshift 

Model Relation 0.0 0.5 1.0 1.5 

Radiative Total 340 190 85 31 

T dy n - A/500 330 186 84 31 

T gas - M 500 332 186 82 31 

Tboi - A/soo 319 151 64 24 

ibol - A/500 317 186 85 31 

£boi-?W 256 95 34 14 

Preheating Total 283 147 59 22 

T dyn - A/500 273 143 56 22 

T gas - A/500 271 143 56 22 

Tbol - A/500 264 134 53 22 

Tbol - A/500 269 143 59 22 

T b oi-T bo i 190 92 48 14 

Feedback Total 342 98 45 13 
Tdyn- A/500 328 96 43 12 
T gas - A/500 327 89 41 11 
Tbol - A/500 305 90 39 10 
Lbol - A/500 339 98 45 13 
L bol -T bol 269 67 32 12 



and ~ 1400 dark matter particles in the (higher resolu- 
tion) Feedback simulation. At z = 0, each model contains 
around 300 clusters above our mass limit, decreasing by 
around an order of magnitude by z = 1.5. 

We also made a number of additional cuts to the cat- 
alogues, specific to each scaling relation. Firstly, we 
noted a small number of systems that were significantly 
offset from the mean relation. On inspection, such ob- 
jects were found to be erroneous as they were subclumps 
falling into neighbouring clusters. Thus, for each rela- 
tion, we discarded all objects with A log(F) > 0.1, larger 
than intrinsic scatter in the Td yn — M500, Tgas — -W500 
and Tboi — M500 relations; and Alog(Y) > 0.5 in the 
Tbol — and Tbol — Tjoi relations, respectively. Sec- 

ondly, for the Tbol — Tbol relation, we made an additional 
cut in temperature, such that the catalogues were com- 
plete in Tboi (excluding those clusters classed as outliers) . 
For the Radiative model, the minimum temperatures are 

fc2hoi,mm = t - 74 ' L0 ' L25 ' L35 l keV; for the Preheatin 9 
model, /cT b oi,mi„ = [0.70,0.96, 1.1, 1.37] keV; and for the 

Feedback model, fcT bo i,min = [0.59,1.12,1.31,1.58] keV, 

for z = [0, 0.5, 1, 1.5]. The numbers of clusters remaining 

in each of the relations after these cuts are also listed in 

Table [2 

4. RESULTS 

4.1. Scaling relations at redshift zero 

We first present the scaling relations at z = as they 
will form the basis for measuring evolution in the cluster 
properties with redshift. The parameters a and Cq(0) 
listed in Table [3 are determined from the best least- 
squares fit to the relation 

logF = logC (0) + alogA, (11) 

where X and Y represent the appropriate data sets 
in units of 10 14 h _1 M 0) 1 keV and 10 42 h~ 2 ergs" 1 for 
mass, temperature and luminosity, respectively. We will 
consider each relation in turn. 

Figure [21 illustrates the Tj yn — M500 relation for each 
of the three simulations at z = 0, with best-fit relations 



Table 2 
Best-fit scaling relations 



Relation 


Model a 


log C (0) logYo 


A 


Tiyn 


- A/500 


Radiative 0.70 


0.34 


0.34 


1.1 






Preheating 0.70 


0.33 


0.33 


1.1 






Feedback 0.69 


0.33 


0.33 


1.2 


Tgas 


- A/ 50 o 


Radiative 0.61 


0.33 


0.33 


0.9 






Preheating 0.61 


0.35 


0.35 


0.9 






Feedback 0.61 


0.35 


0.35 


1.1 


Tool 


- A/ 500 


Radiative 0.59 


0.38 


0.37 


0.5 






Preheating 0.61 


0.35 


0.35 


0.8 






Feedback 0.64 


0.33 


0.33 


1.2 


Tbol 


- A/500 


Radiative 1.82 


1.36 


1.36 


3.9 






Preheating 1.92 


1.40 


1.39 


3.1 






Feedback 2.10 


1.40 


1.40 


3.2 


Z-bol 


— TjoI 


Radiative 3.06 


0.19 


0.20 


1.9 






Preheating 3.05 


0.26 


0.24 


0.7 






Feedback 3.13 


0.28 


0.28 


-0.6 




13.5 14 14.5 

i°g M 5 oo/ h_lM o 



Fig. 2. — The dynamical temperature-mass relation within R500 
at z = 0. Radiative clusters are plotted as circles (with the solid 
line denoting the best-fit relation), Preheating as squares (dashed 
line) and Feedback as triangles (dash-dotted line) respectively. 



overplotted as straight lines. The dynamical temperature 
is dominated by the contribution from the more massive 
dark matter particles, and so the resulting three relations 
are almost identical. The measured slope of the relation 
is a ~ 0.7 (Table|2J, close to, but slightly larger than the 
self-similar value (a — 2/3); this deviation is due to the 
variation of concentration with cluster mass. When the 
mass- weighted temperature of hot gas is used instead, the 
relation becomes flatter than the self-similar prediction, 
with a ~ 0.6. This is expected as the excess entropy 
generation due to cooling and heating is more effective 
in lower mass clusters (MTKP02). 

Shown in Figure0is the Tbol — -M500 relation for each of 
the 3 models. Cool, dense gas dominates Tboi and so this 
temperature is more susceptible to fluctuations caused 
by merging substructure, leading to an increase in the 
scatter when compared to Figure |2 Again, the slope is 
flatter than the self-similar prediction, due to the effects 
of excess entropy. Differences between the models are 
larger than for the dynamical temperature but are less 
than the intrinsic scatter. 

Finally, we consider relations involving the bolometric 
luminosity of the cluster. Fitting the relation between 
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Fig. 3. — The bolometric emission-weighted temperature-mass 
relation within -R500 at 2 = for clusters in the 3 simulations. 
Symbols and lines are as in Figure 121 




log kT b0| /keV 

Fig. 4. — The bolometric luminosity-temperature relation within 
-R500 at z = of clusters in the 3 simulations. Symbols and lines 
are as in Figure l2l 



luminosity and mass, we find a slope in the range a ~ 
1.8 — 2.1, significantly steeper than the self-similar pre- 
diction (a = 4/3). The departure from self-similarity is 
exacerbated when we plot bolometric luminosity against 
temperature (Figure |3J. Here, a ~ 3.1 in all mod- 
els, compared to a — 2 for the self-similar case. The 
Lboi— Tboi relations from the three simulations are in rea- 
sonable agreement with one another and in good agree- 
ment with the observed luminosity-temperature relation 
(see MTKP02,KTJP04). 

In summary, all three models successfully generate ex- 
cess entropy in order to break self-similarity at the level 
required by the observations at low redshift (z ~ 0). 
Thus, based on the local scaling relations alone, we can- 
not easily discriminate between the source of the entropy 
excess in clusters: whether it is mainly due to radia- 
tive cooling, additional uniform heating at high redshift 
(prior to cluster formation) or localized heating from 
galaxy formation at all redshifts. 

4.2. Evolution of scaling relations with redshift 

We now examine whether this degeneracy between 
models in the scaling relations at z = can be broken 
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Fig. 5. — The normalization of the various temperature-mass re- 
lations as a function of redshift, for clusters in the Radiative (solid 
line), Preheating (dashed line) and Feedback (dot-dashed line) 
simulations: dynamical temperature (top panel), mass- weighted 
gas temperature (middle panel) and bolometric, emission-weighted 
temperature (bottom panel). 



by examining the cluster population at higher redshifts 
(z = 0.5, 1, 1.5). None of the relations require a signif- 
icant variation in a with redshift. To make our results 
easier to interpret, therefore, we use simple power-law 
relations of the form given in equation 0] with a fixed at 
the z — values given in Table [3 

To find the evolution of each relation, we first deter- 
mine the normalizations, Co, and their corresponding er- 
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Fig. 6. — The normalization of the bolometric luminosity-mass 
relation as a function of redshift, for clusters in the Radiative (solid 
line), Preheating (dashed) and Feedback (dot-dashed) simulations. 



ror bars, at each redshift in the same manner as described 
for redshift zero in Section 4.1 above. We then minimize 
the x~ squared to obtain parameters Y$ and A as listed 
in Table to fit the relation 

logC = lo g r +Alog(l + z). (12) 

4.2.1. Temperature-Mass Evolution 

In Figure^] we present values of log(C ) versus redshift 
for the three temperature-mass relations, with the best- 
fit straight line ovcrplotted. For the Td yn — -M500 relation 
(upper panel), we find similar evolution parameters for 
the three models, A = 1.1 — 1.2, confirming that including 
the effects of baryonic physics does not significantly affect 
cluster dynamics. The slight excess over the self-similar 
value of A = 1 is consistent with the changing cluster 
concentrations. 

However, both the mass-weighted temperature (mid- 
dle panel) and especially the emission-weighted temper- 
ature (lower panel) show significant variation between 
the three models. In each case the Feedback simulation 
approximately follows the scaling found for the dynam- 
ical temperature, with the Preheating and the Radiative 
simulations showing progressively larger deviations be- 
low the expected normalization as the redshift increases. 
The explanation for this lies in the variation of temper- 
ature of gas particles within each cluster and how this 
changes with redshift in the different models. We shall 
explore this further in Section [SJ 

4.2.2. Luminosity- Mass Evolution 

Figure[f)lillustrates the normalization of the Lboi— M500 
relation versus redshift for all three models. The Pre- 
heating and Feedback models evolve almost identically 
with redshift (A ~ 3), but the Radiative run evolves 
more strongly (A ~ 4). These bracket the self-similar 
value, A — 3.5, however, this agreement is somewhat 
coincidental because the slope of the relation at fixed 
redshift is much steeper than expected (a ~ 1.8-2.1 
rather than 1.3). The reason why the Radiative simu- 
lation has steeper evolution is because of enhanced emis- 
sion from cool gas at high redshift relative that that at 
low redshift — see discussion in Section [5] 
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Fig. 7. — The bolometric luminosity— temperature relation for 
clusters at z = 0,0.5,1,1.5 (solid, dashed, dot-dashed, dotted lines 
respectively), various redshifts. The top panel is for the Radia- 
tive simulation, the middle for the Preheating simulation and the 
bottom for the Feedback simulation. 



4.2.3. Luminosity- Temperature Evolution 

Finally, we consider the evolution of the Lboi — Tboi 
relation, with the relations at each redshift shown ex- 
plicitly for each model in Figured and the variation of 
normalization with redshift illustrated in Figure [S] It is 
interesting to note that the values of A are significantly 
different between all three models: the Feedback model 
predicts mildly negative evolution (A = —0.6), the Pre- 
heating mildly positive evolution (^4 = 0.7) and the Ra- 
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log (1+z) 

Fig. 8. — The normalization of the bolometric luminosity- 
temperature relation as a function of redshift, for clusters in the Ra- 
diative (solid line), Preheating (dashed) and Feedback (dot-dashed) 
simulations. 



diative strongly positive evolution (A — 1.9). The latter 
two models straddle the self-similar value (A — 1.5). 

The difference in slopes between the Feedback and Pre- 
heating runs is driven by the differences in their tempera- 
ture. The further difference between the Preheating and 
Radiative runs comes roughly equally from the tempera- 
ture and luminosity evolution. 

5. DISCUSSION 

In this paper, we have focused on the evolution of clus- 
ter scaling relations in three simulations, each adopting a 
different model for non-gravitational processes that affect 
the intracluster gas. In the first, Radiative model, the gas 
could cool radiatively and a significant fraction cooled 
down to low temperatures and formed stars (MTKP02). 
In the Preheating model, the same was true, although the 
gas was additionally heated uniformly and impulsively by 
1.5 keV per particle at z = 4, before cluster formation. 
In the third, Feedback model, the heating rate was local 
and quasi-continuous, in proportion to the star-formation 
rate from cooled gas. All three models are able to gen- 
erate the required level of excess core entropy in order 
to reproduce the L x - T x relation at z = (MTKP02, 
KTJP04). 

The most striking result presented in this paper is that 
the three models predict widely different Lx — Tx rela- 
tions at high redshifts. The Radiative model predicts 
strongly positive evolution, the Preheating model mildly 
positive evolution and the Feedback model, mildly nega- 
tive evolution. At this point, it should be stressed that 
the values of A presented in Table |3 should not be taken 
too seriously. No attempt has been made to convert 
the bolometric, emission-weighted fluxes used in this pa- 
per to observable X-ray fluxes in different instrumental 
bands. Also, the volume of our simulation boxes is such 
that we only get a modest number of relatively poor 
clusters at high redshift. Nevertheless, the qualitative 
difference between the models is very encouraging and 
suggests that evolution of X-ray properties may act as a 
strong discriminant between mode ls in the futu re. 

.Previous work (e.g. iPearce et all 12005 

iMuanwong et al. 1 I200H) has made great play of the 



fact that radiative cooling can remove low-entropy ma- 
terial and lead to a raising of the gas temperature above 
the virial temperature of the host halo. That effect is 
reproduced by the Radiative simulation in this paper, 
but it is interesting to note that the bolometric X-ray 
temperature exceeds the dynamical temperature of the 
clusters only at very low redshift, z < 0.1. At higher 
redshifts it falls below the dynamical temperature and 
is a factor of 1.6 lower by z — 1.5. This departure from 
self-similarity is a consequence of the changing density 
parameter, Q, in the concordance ACDM cosmology: 
at high redshift £1 is close to unity and structures grow 
freely; at lower redshifts f2 falls well below unity and the 
rate of growth of cosmic structures declines. 

The behavior of the Preheating simulation is similar, 
although the relative decline in the ratio of the gas 
to the dynamical temperature is smaller. The effect 
cannot, therefore, be due to the cooling of intraclus- 
ter gas in the cluster cores between a redshift of 1.5 
and the present — in the Preheating run very little gas 
cools below a redshift of 4 and so that cool core gas 
would still be there today. Instead, we attribute the 
presence of cool gas to the accretion of low temper- 
ature subclumps. Suc h accretion is a ubiquitous fea - 
ture of clusters (e.g. iRowlev. Thomas fc Kav I [20041 . 
Cool gas is seen in maps of clusters at low redshift 
l)Onuora. Kav fc Thomas ll200.1l) and would be expected 
to be much more prevalent in clusters at high redshift in 
the ACDM cosmology. 

To test this hypothesis, we measured the temperature 
variation within each cluster as follows. First we aver- 
aged properties within 20 spherical annuli out to -R500, to 
create smoothed dynamical, Td yn , and gas temperature, 
r gas , profiles. Then we measured the mean deviation (in 
log space) of the gas temperature from the local dynam- 
ical temperature 

T = ( lo Sio T * ~ l °Sw T dyn) , (13) 

i 

and the root-mean square deviation of the temperature, 
(Jt-i from the mean, 

°r = Jj E ( lQ gio T i ~ lo Sio r gas ) 2 , (14) 

i 

where N = and the sum runs over all hot gas particles 
(Ti > 10 5 K) within #500. 

As an example, FigureElshows the values of r (crosses) 
and cjt (circles) for each cluster at z = 1 in the Radia- 
tive simulation. This particular example has been chosen 
simply because the two properties are well-separated and 
easy to distinguish on the plot. As can be seen, the mean 
gas temperature of the more massive clusters typically 
lies below the local dynamical temperature, by as much 
as 15 per cent at this redshift. However, the dispersion 
in temperature is much larger, typically a factor of 1.5, 
so there will be fluctuations both above and below the 
dynamical temperature. 

Figure ITU1 demonstrates visually the evolution of r and 
<7t with redshift. At each redshift, the plot shoes the av- 
erage value of r over all clusters with masses greater than 
1.2xl0 13 /i -1 Mq. The half-width of the shaded regions 
represent the average values of ctt, divided by 10 for clar- 
ity. Concentrating first on the Radiative simulation, it 
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Fig. 9. — Values of r (crosses) and ctt (circles) for individual 
clusters in the Radiative simulation at z = 1. 
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Fig. 10. — The mean value of r for clusters in the Radiative (solid 
line), Preheating (dashed line) and Feedback (dash-dotted line) sim- 
ulations as a function of redshift. The half-width of the shaded 
regions corresponds to the mean value of ctt at that redshift. All 
averages are for clusters with mass greater than 1.2x10 h -1 Mq. 



can be seen that the mean cluster temperature increases 
relative to the virial temperature over time and that the 
dispersion in temperature decreases. This is consistent 
with a decreasing amount of substructure within the clus- 
ters at lower redshift, although it should be noted that 
part of the effect is due to the narrower range of cluster 
masses resolved by the simulations at high redshift, as 
the average value of r decreases with increasing cluster 
mass. 

The behavior of the Preheating simulation mimics that 
of the Radiative one, but with a bias to higher mean tem- 
peratures. The Feedback simulation, however, is quite 
different. It shows a much larger dispersion than the 
other two, but no bias to low temperatures at high red- 
shift. That is because gas is free to cool down to low 
temperatures but some of that gas is then heated back up 
high temperatures by the feedback. At low redshift cool- 
ing becomes less important and the Feedback run then 
shows a slight rise in r, matching that seen in the other 
two runs. 

Our results have two very important implications for 
observations of high redshift clusters. Firstly, the behav- 



ior of the Lx — Tx relation at high redshift will determine 
the number of high redshift clusters to be found in sur- 
veys such as the on going XMM-Newton Cluster Survey 
(R omer et al. 11200 if) and this will have a significant im- 
pact upon their use in probing cosmological parameters. 
A positive evolution such as that shown by the Radia- 
tive simulation will yield many more observable high- 
rcdshift clusters than the negative evolution of the Feed- 
back model. As discussed in Section [21 and summarized 
in Figure ^ the observational situation is far from clear 
but does seem to indicate positive evolution. 

Turning this argument around, our results suggest that 
observational constraints on the degree of evolution of 
the Lx — Tx relation will allow interesting constraints to 
be placed on the source of entropy generation in clus- 
ters, in particular the relative role of cooling and heating 
and whether most of the heating of the intracluster gas 
occurred at high redshift (as in the Preheating model) 
or was a continuous function of redshift (as in the Feed- 
back model). Taking our results at face value with recent 
observations would suggest that our Feedback model is 
generating too much excess entropy at z < 1.5 and that 
the bulk of the heating must have occurred at higher red- 
shift. However, we stress once again that this result is 
very tentative. 

6. CONCLUSIONS 

The evolution of X-ray cluster scaling relations are a 
crucial component when constraining cosmological pa- 
rameters with clusters. Observational studies at low red- 
shift have already shown that the scaling relations de- 
viate from self-similar expectations, attributed to non- 
gravitational heating and cooling processes, but their 
redshift dependence is only starting to be explored. In 
this paper we have investigated the sensitivity of the X- 
ray scaling relations to the nature of heating processes, 
using three numerical simulations of the ACDM cosmol- 
ogy with different heating models. While all three simu- 
lations reproduce more or less the same scaling relations 
at z — (as they were designed to produce the correct 
level of excess entropy), they predict significantly dif- 
ferent results for the evolution of the Lx-Tx relation to 
z = 1.5. 

In conclusion, our findings strongly suggest that the 
relative abundance of high and low redshift clusters 
will place interesting constraints on the nature of non- 
gravitational entropy generation in clusters. First indi- 
cations are that an early and widespread preheating of 
the ICM is to be preferred to an extended period of pre- 
heating that is associated with galaxy formation. How- 
ever much more detailed modeling is required and the 
observational picture is as yet unclear. 
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